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Abstract 

A recent theory that determines the properties of disordered soHds as the sohd 
accumulates damage is apphed to the special case of fiber bundles with global load 
sharing and is shown to be exact in this case. The theory postulates that the prob- 
ability of observing a given emergent damage state is obtained by maximizing the 
emergent entropy as defined by Shannon subject to energetic constraints. This theory 
yields the known exact results for the fiber-bundle model with global load sharing 
and holds for any quenched-disorder distribution. It further defines how the entropy 
evolves as a function of stress, and shows definitively how the concepts of temper- 
ature and entropy emerge in a problem where all statistics derive from the initial 
quenched disorder. A previously unnoticed phase transition is shown to exist as the 
entropy goes through a maximum. In general, this entropy-maximum transition oc- 
curs at a different point in strain history than the stress-maximum transition with 
the precise location depending entirely on the quenched-disorder distribution. 
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1 INTRODUCTION 



The fiber-bundle model with global load sharing is a simple model for failure in 
tension introduced almost 80 years ago [1-3] and having received considerable 
attention and extensions over the past 15 years [4-19] . Although this model 
may have little pertinence to real fibrous systems such as a rope breaking in 
tension, it is of interest because it possesses exact analytical properties. 
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In recent work [20-22] , we have developed a general statistical theory for de- 
termining the properties of a disordered solid that is accumulating irreversible 
damage due to cracking under stress. The ensembles in this theory are cre- 
ated by considering different realizations of the quenched disorder in the local 
breaking strength of the material. Such realizations are made either for the 
system as a whole, or of more pertinence to real systems, by dividing a given 
system into smaller "mesovolumes" and letting each mesovolume correspond 
to a different realization of the quenched disorder. The ensembles so obtained 
have nothing to do with thermal fluctuations (molecular dynamics). Given a 
system with such quenched disorder, our theory determines the probability of 
emergent crack states by maximizing Shannon's measure of disorder subject 
to constraints coming from the energetics of the fracture process. 

In the present paper, we apply this theory to the specific problem of fiber bun- 
dles with global load sharing and demonstrate that for any quenched-disorder 
distribution, it yields the known exact results. Furthermore, a previously un- 
noticed phase transition is demonstrated to exist where entropy goes through 
a maximum. This phase transition is distinct from the well-known stress maxi- 
mum transition and was noticed in the present theory because of the prominent 
role played by entropy. 

However, the importance of our theory is not that it yields a new result in 
this old model, but that it applies and yields analytical results about phase 
transitions for any quasi-static damage model; albeit, for models involving 
crack interactions, approximate treatment of the functional integrations may 
be required to obtain analytical results (such as renormalization or mean-field 
approximations). Using our theory, we recently treated the problem of how the 
mechanical properties of rocks change due to cracks arriving and interacting in 
compressive shear [20-22]. We analytically demonstrated that the localization 
transition observed in experiments is a second-order phase transition when it 
occurs. 



2 THERMODYNAMICS 



Our theory was originally developed for a disordered solid under the influence 
of a stress tensor. For fiber bundles, a much simpler scalar theory applies and 
so in this section, the formalism is rederived in this simplified context. 

A fiber bundle is depicted in Fig. 1. A collection of fibers are stretched 
between two rigid supports. One support is held fixed, while the other is free 
to displace. A load is applied to the bundle through the free support so 
that the fibers are in a state of tension. In this paper, the load F/v will always 
be normalized by N to define the overall tension r = F^/N. Each non-broken 
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Fig. 1. A bundle of N fibers stretched between two rigid supports with a load -Pjv 
applied through the free support. 

fiber in a bundle has the same length L. If, when t — 0, this length is Lg, then 
the measure of strain is s = L/ Lg — 1. Experiments may be performed on the 
bundle either by controlling r or e. 

Each of the fibers has the same Young's modulus which is taken to be unity so 
that the axial strain e of each fiber is identical to the tension in the fiber. The 
fibers have strengths ei, 62,... which are independent random variables 
sampled from a distribution p{e), whose cumulative distribution is defined 
P{e) = jQp{x)dx. As the strain e of the bundle is increased, the individual 
fibers will break once their tension (strain) gets to their fixed strength thresh- 
old. The assumption of "global load sharing" is that when one of the fibers 
breaks at fixed load, all the other fibers will extend by the same amount thus 
increasing the tension in each of the surviving fibers so that the load as a 
whole is always supported entirely by the surviving fibers. All of this defines 
the fiber bundle model with global load sharing. Of interest are the mechani- 
cal properties of such bundles as averaged over all possible realizations of the 
fiber strengths. 

We first need to know the probability pj of observing one of the realizations 
to be in a particular state j of damage when the ensemble as a whole is at 
an applied strain s. In the fiber bundle model, a damage state j is defined 
by which of the fibers are broken. One could define j using a local order 
parameter that is 1 if a fiber is intact and if the fiber is broken. 

Our theory postulates that the fraction pj of all realizations observed to be in 
state j is obtained by maximizing Shannon's measure of disorder 

S^~j:P^lnp, (1) 
j 

subject to constraints. Such constraints must involve the independent variables 
of S. To identify the independent variables, we consider how both S and the 
average energy in the ensemble of bundles change as the strain is increased. 

When e increases to e + de, there is both a work carried out in reversibly 
stretching the fibers and an additional work carried out due to irreversible 



fiber breaks. Due to breaking, some of the members of tlie ensemble (individual 
realizations of the disorder) will be led out of their current damage state and 
into state j, while others that were in state j will transfer to still different 
states. If there is a difference in the number of members entering and leaving 
state j, there will be a change dpj in the occupation probability of state j and 
such changes are what cause Shannon's disorder measure S to change. 

The average energy density (average energy normalized by N) in the ensemble 
is given hy U — J2jPjEj. Here, Ej is the energy density required to create state 
j at imposed strain e and as averaged over all members that have been led 
to state j. Depending on the breaking strengths of a given realization, the 
work performed in arriving at state j can be different. It is through Ej that 
all dependence on the quenched-disorder distribution enters the problem. The 
change that occurs when e increases to e + de is 

dU^J2Ejdp^ + J2PjdEj. (2) 

j j 

The first term is the energy expended in changing the disorder over the col- 
lection of realizations. It is thus proportional to the disorder change and can 
be written 

TdS = ^Ejdpj. (3) 
j 

The second term is written 

fde^Y.P3dEy (4) 

3 

and represents both the reversible stretching energy in those members that 
did not experience breaks during the deformation increment, as well as the 
irreversible energy changes due to all the breaks that did not result in a net 
change in the occupation numbers of each state. 

This decomposition of the energy increment can be thought of as follows. 
When breaks occur throughout the ensemble of realizations as e increases to 
e-\-de, there is a fiow of members between the states. This fiow involves energy 
changes due to fibers breaking in the interval de. It may be resolved into a 
uniform "incompressible" part whose associated energy is included within / de 
as well as a non-uniform "compressible" part associated with more members 
arriving in a given state than are leaving that state. The energy associated 
with the non-uniform fiow between states is entirely contained in TdS. 

From these expressions it may be concluded that if U is to be treated as a 
fundamental function, then U = U{S,e), or equivalently if S is to be treated as 
the fundamental function then S = S{U,e). In other words, the independent 
variables that must be involved in the constraints on the maximization of S 



4 



are U and e. The proportionality constants T and / are defined 



The state function / is something different than the overall tension r since we 
also have that t de — dU . Thus, in general, {t — f)de — TdS so that f ^ t due 
to fibers breaking in a positive increment de. If strain were to be decreased, 
fibers do not break and so dS = Q and the state function / would be defined 
using only the purely elastic part of the energy changes dEj. Changes in / in 
this case are equivalent to changes in r. Last, since we have taken S to be 
extensive (proportional to N) while [/ is a density independent of N ^ T comes 
out being proportional to 1/N . This choice is made so that factors of N do 
not clutter the equations that follow. 

The constraint involving e is that each non-broken fiber throughout the entire 
ensemble has the same length which implies Sj — e. The constraint involving 
U is that U — ^jPjEj. Carrying out the maximization of S subject to these 
constraints using Lagrangian multipliers gives the probability distribution as 

where /? = 1/T and where the partition function Z is defined 

Z{P,e)^^e-f'^^^^\ (7) 
j 

The free energy F = F{P, e) in this ensemble (for this set of constraints) is the 
Legendre transform of U with respect to S; i.e., F = U — TS. Differentiation 
then gives 

dF ^ p-'^Sdp + fde. (8) 

Upon introducing U — J2jPjEj and S — — ^jPjlnpj into the Legendre trans- 
form F ~U — S/ (5 and using that Inpj = —^Ej — In Z we obtain the standard 
result 

pF = -\n Z. (9) 

Thus, the standard canonical ensemble emerges in this problem where quenched 
disorder alone (not molecular fiuctuations) is responsible for the existence of 
ensembles. 

An important question in such an approach is whether anything precise can be 
said about the temperature T = Indeed, (3 can formally be found as the 
solution of a differential equation. This differential equation is obtained from 
the previously stated but unused fact that dU — rde which can be written 

dU X - dpi ^ . ^ dE^ , , 

-^=i:f^.+E«^. (10) 

3 3 
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Fig. 2. A strain-controlled experiment. The solid line is the load path followed during 
the experiment. The dashed lines represent the path that would be followed if strain 
were to be decreased at some point during the experiment. When a fiber breaks at 
constant s, the load r must be reduced as represented by the vertical drops. 



Now from Eq. (6), we obtain that 



dPj 
de 

while from Eqs. (8) and (9) 
dlnZ 



dEj ^ _ ^ dp _ dlnZ 
ds ^ de de 



(11) 



de 



= -F 



d/3 
de 



de 



-U 



dp 
de 



Pf. 



(12) 



Since each member is at the same e, each member has its own tj, and so 
T = J2jPjTj- We then obtain the differential equation for /3 in the form 



d3 

a-r- + 6/3 + c = 
de 



with coefficients given by 



^p,E,{U-E,); b = J2PjEjlf 



dEj^ 
de 



^ (dE. 



(13) 



(14) 



Since pj = pj{P,e), this equation is non-linear and thus difficult to solve. In 
the present work, we demonstrate that a proposed function P — P{e) exactly 
satisfies this equation and thus is the true "fiber-bundle" p. To make progress, 
we next obtain Ej and Tj for the specific problem of fiber bundles with global- 
load sharing. 



3 FIBER BUNDLE MODEL 



Figure 2 details the history of how the overall tension r might evolve when con- 
trolled variations in e are applied to a bundle having fiber strengths £1, £2---£jv- 
In this particular example, n = 3 fibers have broken when the strain is at e. 
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The load on the bundle is equally shared by the N — n surviving fibers so that 



r„ = ^(l-^) (15) 

which is a relation independent of the history; i.e. it depends only on the 
actual state of the bundle through the number of broken fibers n, and not on 
the breaking thresholds £i, £2 

The total work density Ej for states j consisting of n broken fibers is the area 
under the charging curve in a particular realization (the area under the solid 
line in Fig. 2) 

where from Eq. (15) we have used t^(x) = x{l — m/N) and where by con- 
vention £„+i = £ is the final applied strain. A direct recursion gives exactly 

n\ -J^ el 




The first term here is the elastic energy that can be reversibly recovered upon 
decreasing the strain while the second term represents the energy irreversibly 
consumed in the breaking process. Both contributions can be directly vizual- 
ized in Fig. 2. 

Equation (16) is next averaged over the quenched disorder to obtain the av- 
erage energy Ej needed to create state j. Each breaking threshold em is an 
independant variable, randomly distributed according to p{ern) under the con- 
dition that < < They are therefore distributed according to the prob- 
ability density p{em)/P{e) where the normalization factor accounts for the 
fact that the upper limit e is independent of the threshold values e^,- Thus, 
averaging over the quenched disorder gives 

1 /"^ x' 




r %-p{x)dx (17) 

Jn 1 



P{e)Jo 2 

where h{e) designates the average energy that is lost when each fiber breaks. 

The average work density (Hamiltonian) required to create the state j, aver- 
aged over all realizations of the quenched disorder is then in general 



(18) 



For example, under the special assumption that the breaking strengths are 
randomly sampled from a uniform distribution on < ^ 1, we have that 
p{e) = 1, P{e) = e, h{e) = e^G, and Ej = 6^2 - {rij/Ny/S. 
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4 Average Properties 



We now apply the above theory and determine the thermodynamic functions 
as a function of applied strain e. The following analysis is valid for any properly 
normalized quenched-disorder distribution p{s). 

The cumulative distribution P{s) = Jq p{x) dx is the probability that any one 
fiber has broken when the strain is at e. Thus, the fraction of all possible 
realizations having rij broken fibers and N — rij unbroken fibers is exactly 

^exact = (^l — p^^-^j p"-j_ ^19^ 

This distribution may be written 



^exact ^ 



—rij In 



(20) 



where po = (1 — -P)^ is the probability of the unbroken state j = 0. Upon using 
the Hamiltonian of Eq. (18), the postulate of entropy maximization predicts 
this same distribution to be given by 



Pi =Poexp 




(21) 



where po — exp{—(3e^/2)/Z is the probability of the unbroken state. 

These two distributions are both Boltzmannians in the number nj c 
fibers and are identical if the temperature T — 1/P is given by 




ey2-h{e 



(22) 



If it can be shown that this (3 satisfies the differential equation of Eq. (13), then 
our theory is exact when applied to fiber bundles with global load sharing. 

From the definition P{e)h{6) = Jq p{x)x^ /2 dx, one has h{e) < for any 
distribution p(£). Thus, /3(e) is a negative increasing function up to the strain 
point £ — £/3 where it smoothly goes to zero. The inflection point is obtained 
from the condition that P(£fl) = 1/2 and defines a previously unnoticed phase 
transition that will be shown to be distinct from the transition at peak stress. 
When e > e^, (3 becomes a positive increasing function of £. 

There are two key average properties upon which all the thermodynamic func- 
tions depend; namely, the average fraction of broken fibers in each bundle 
{rij/N) and the average of this fraction squared {{rij /NY). Using the exact 
probabilities of Eq. (20) [which is equivalent to using the /5 of Eq. (22) in our 
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probability law], one obtains 

I) = Ep.| = E -n"^^"(l - ^')"-" (23) 

where — A^!/[n!(A^ — n)!] defines the number of ways of selecting n objects 
from a group of distinguishable items. The binomial theorem states that 

N 

(^ + Z/)'' = Ec^^V-"- (24) 

n=0 

Upon differentiating this equation with respect to x and then multiplying by 
x/N, gives that when x = P and y = 1 — P 

'fl = P (26) 

which is a known result consistent with the meaning of P. Differentiating a 
second time yields 

Using these two results, the other averages defining the thermodynamic vari- 
ables are easily read off. 

The average stress T{e) is thus obtained to be 

^ = HVj^j = (1 - (27) 



which is initially an increasing function of e having the slope 

%^1-P-ep. (28) 
as 

This slope goes to zero and the stress has a maximum at any strain point £r 
that is a solution of 1 — P{£t) — £tP{£t) = 0; i.e., at the point(s) where 

ErPiSr) = 1 - P(x) dx (29) 

Jo 

admits a solution. This is a known exact result [4,10]. In general, we can 
expect that Sr^ep. The condition required for equality of these strain points 
is that simultaneously £p{e) = 1/2 and JqP{x) dx — 1/2 which for monotonic 
distribution functions p{e) can only occur with the uniform distribution p{e) — 
1 in which case Sr = = 1/2 and the two transitions coincide. If p{e) is a 
monotonic increasing function of £, then Sr < £f3 while if it is a decreasing 
function of e, then Et > ep. For non-monotonic distributions, there can be an 
arbitrary number of stress maximas and Eq. (29) can have either no solutions 
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or multiple solutions [10] . The nature of the phase transitions at the distinct 
strain points Ejs and is discussed in the following section. 



The average energy in the ensemble is 



U^J2pjEj^{l-pf^ + hP. (30) 
j ^ 

Again recalling the definition of h from Eq. (17) gives 

dlJ , d(hP) , , 

which is also the equation that gives the differential equation for temperature. 
This is sufficient for demonstrating that the 13(e) of Eq. (22) satisfies its differ- 
ential equation. Nonetheless, as a consistency test, the coefficients a, b, and c 
defined in Eq. (14) will be derived and the differential equation will explicitly 
be shown to be satisfied. 

To obtain the state function /, we need first 



dEi Tin 




h 



de N 

From Eq. (4) and the lemma of Eq. (25) we then have 



(32) 



f-EPj^ = i^-P)s+pi"-^-h]. (33) 



, de ^ ^ 



The variation of the entropy is obtained from the energy balance as 

Together with the initial condition S'(O) = 0, this is readily integrated to give 

5= -iV[PlnP + (l-P)ln(l-P)]. (35) 

This expression is the classical Shannon result for a set of N random variables 
in two possible states having probabilities P and 1 — P, which is precisely 
the case of the fiber bundle with global sharing. This is another consistency 
check. This entropy S increases from zero [total certainty that every member 
is intact] and goes smoothly through a maximum [total uncertainty as to what 
state a member may be in] at the same strain point P{ep) — 1/2 where 13 goes 
to zero. After the smooth maximum, S decreases to reach zero if P{s) reaches 
1 in which case there is total certainty that each member is entirely broken. 
We note that Eq. (35) can also be directly obtained from the Shannon formula 
upon applying the binomial theorem. 



10 



Finally, the free energy F is again obtained from its Legendre transform defi- 
nition F = U - S/p to be 



\ In(l-P) 
^2 y In(l-P) -InP' 



(36) 



At the point e/3 defined by P^Sp) = 1/2, the free energy diverges due to the 
fact that I3{efi) — while S{e^ remains finite. So long as 7^ e^, the free 
energy does not diverge when (if) e — 8^-. In passing, we also note that 



(1 - Pfp-^'l'' 



N/{e^/2-h) 



(37) 



is the exact expression of the partition function. 

With the above results established, we now obtain the coefficients a, 6, c of the 
differential equation for the temperature as 



N \ 



l^(ph^-'l\(pe-/l+ph); 



(38) 



Using these together with Eq. (22) for (5 and its derivative 



dl3 
de 



Np 



P(l - P)(h - ey2) 



+ N\n 



1-P 
P 



e + {h - sy2)p/P' 



(/l- £2/2)2 



(39) 



shows that the differential equation a dp / de + b(3 + c — is exactly satisfied. 



5 Phase Transitions 



5.1 The entropy maximum transition 



The strain point Sf^ defined by P^Sf)) — 1/2 is where simultaneously (3 = 0, 
the entropy is a maximum, and the free energy diverges. It is distinct from 
the stress-maximum transition(s) at Er- The interpretation of £ = 6/3 as a 
phase transition is natural, since the most probable configuration of a bundle 
suddenly jumps from being entirely intact to entirely broken. What are the 
measurable manifestations of this transition at 

Since the entropy is a maximum at this point, the fiuctuations between realiza- 
tions should also be a maximum as we now demonstrate. Define the fraction of 
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broken fibers in state j to be pj = n.j/N, and define the average of this fraction 
to be p = {pj) = P where P = P{e) is again the probability of a fiber being 
broken. The quantity of interest here is the root-mean-square fluctuation Ap 
in the fraction of broken fibers given by 



Ap=y(Z^=^(pj)-p^=VfJ^ (40) 

which indeed goes through a maximum at = 1/2 as expected. This max- 

imum is something that can be directly measured in numerical experiments 
on fiber bundles but has never before been commented on. The reason it has 
been discovered in the present theory is because entropy and temperature are 
explicitly present. We encourage someone to numerically measure Ap and to 
verify that it is a maximum at the transition point e/3. Recall that for mono- 
tonic quenched-disorder distributions, if p{e) is a decreasing function (more 
weak fibers than strong fibers), then < Et-. So the numerical observation 
in this case should be that Ap goes through a maximum prior to peak stress. 
Equivalent comments hold when p{e) is an increasing function and S/^ > Et-. 

The other fluctuations that are potentially of interest include the root-mean- 
square stress fluctuation 

Ar=^[(A^=^{rj)-T- = eAp (41) 
and the root-mean-square energy fluctuation 

AU = ^/{A^ = ^{Ef)-U- =i$-^) ^^^^ 

but since these are simply proportional to Ap it seems that the interesting 
signature of this phase transition is the maximum in Ap. 

The "order" of this transition is not classically deflned in the Ehrenfest scheme; 
however, it seems inappropriate to call it a continuous transition because the 
free energy is singular at Sp as shown above. Nonetheless, in the limit as £ ^ £^ 
the singular part of F diverges according to the scaling law 

\n2 [el-h{ep)/2] 
^'-—8 ^) ^'-''^ ^^^^ 

which has a (trivial) universal exponent. Since F and its derivatives are diffi- 
cult to numerically measure, the principle manifestation of this phase transi- 
tion remains Ap going through a maximum at a point Sr- 
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5.2 The stress maximum transition 



The phase transition at peak stress £ = is the one that researchers up to now 
have focused on. Since S and its derivatives remain continuous and finite there, 
an attempt to classify this transition as being first-order or second-order is 
meaningless. The label of "continuous" transition seems the most appropriate. 

There is universality at this transition due to the fact that for any monotonic 
analytic quenched-disorder distribution. t{£) is a smooth analytic function so 
that upon developing this function in the neighborhood of its maximum one is 
guaranteed |r — t(£t-)| ~ |e — St-P with the exponent of course being indepen- 
dent of the distribution p{s) or other model parameters. Such development of 
an analytic function about a critical point is the way any mean-field theory 
acquires a universal scaling law. However, Kloster, Hansen and Hemmer [10] 
do demonstrate that a non-analytic quenched-disorder distribution can lead to 
highly non-analytic stress-strain behavior that thus falls outside the quadratic 
universality class. Also, in other damage models in which elastic interaction 
between cracks (damage points) is important, the simple observation of a 
model exhibiting an averaged stress maximum does not by itself guarantee an 
exactly quadratic stress-strain relation in the neighborhood of the maximum. 
Such a relation may be non-analytic due to a diverging correlation length in 
the correlation between cracks. 

The quadratic stress maximum means that de/dr diverges as |r — T{er)\~^^'^ 
[11,19] . Accordingly, the average rate at which the fraction of broken fibers 
increases with stress dp/dr = pde/dr also diverges as |r— r(£^)|~^/^ [4,15] . Us- 
ing additional considerations not developed in this paper, one can futher show 
that the average size of avalanches also diverges as |t — T(£r)|~^^^ [10,15,17], 
and that the distribution D{A) of the size A of the avalanches scales as 
D(A) ~ A^^/^ [6,7,10,17] at the stress maximum. These are the principal 
observable characteristics of the stress-maximum transition. 



6 Summary and Conclusions 

Two principal results have been obtained in this paper. First, it was demon- 
strated that for any quenched-disorder distribution used in the fiber-bundle 
model with global load sharing, the probability of the emergent damage states 
may be exactly calculated by maximizing Shannon's entropy subject to con- 
straints. It is through the constraints that the nature of the quenched-disorder 
distribution enters the macroscopic thermodynamics. 

Second, a previously unnoticed phase transition occurs in the fiber-bundle 
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model with global load sharing when the entropy goes through a maximum. 
This phase transition is distinct from the stress-maximum transition; although, 
for the uniform quenched-disorder distribution, the two transitions coincide. 
The principle manifestation of this transition is that the root-mean-square 
fluctuation in the number of broken fibers will go through a maximum, which 
is a quantity that can be directly measured in numerical experiments. 

To conclude, we postulate that the probability of emergent states can always 
be calculated through entropy maximization for any damage model of interest 
including models in which there is elastic interaction between the local dam- 
age (order) parameters. Damage models with order-parameter interactions 
are not normally considered amenable to analytical treament; however, using 
our approach a rather standard canonical-ensemble partition function emerges 
and the various functional integration procedures available for studying the 
partition function in the neighborhood of critical points may be employed. 

Such generality is the principal utility of our approach. A concern is the ability 
to produce an expression for the model temperature. The temperature can be 
found in principle as the solution of a well-posed initial-value problem. Un- 
fortunately, the differential problem is highly non-linear and thus difficult to 
solve. Because stress maxima and entropy maxima do not normally coincide, 
a stress-maximum transition may be studied by simply assuming the temper- 
ature to be well behaved in the neighborhood of the stress maxima. But if an 
explicit expression for the temperature is desired, the following approach can 
be employed. 

In models involving crack interactions, there are always a subset of dilute states 
in which the interactions are negligible. The exact probabilities of such states 
can usually be determined from the quenched-disorder distribution alone. A 
comparison between the probabilities calculated from entropy maximization 
and such exact probabilities then determines the model temperature to be 
used in the Boltzmannian. With the temperature so defined, and with a proper 
model Hamiltonian in hand, the partition function can be analyzed. 
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